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Abstract 

In supersymmetric theories, the top quark can decay into its scalar partner plus a 
neutralino, with an appreciable rate. We calculate the 0(a s ) QCD corrections to this 
decay mode in the minimal supersymmetric extension of the Standard Model. These 
corrections can be either positive or negative and increase logarithmically with the 
gluino mass. For gluino masses below 1 TeV, they are at most of the order of ten 
percent and therefore, well under control. 
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1. Introduction 



It is well-known that in supersymmetric theories (SUSY), the heavy top quark can decay at 
tree-level into its lightest scalar partner and a light neutralino state [fll This possibility 
is intimately related to the large value of the top mass mt- the mixing between the left— 
and right-handed scalar partners of the top quark, and £r, is proportional to m t ; after 
diagonalization of the mass matrix, the lightest stop mass eigenstate t can be much lighter 
than the top quark and all the scalar partners of the light quarks j0J. 

The experimental bounds [[| on the lightest stop and neutralino masses leave a large room 
for this decay to occur. Moreover, light stop squarks and chargino states [and hence the light 
neutralino Xi which is the lightest supersymmetric particle] might explain the large deviation 
of the experimentally measured value of the ratio Rb = T(Z — > bb)/T(Z —>■ hadrons) from 
the Standard Model prediction M. In fact, in the area of the MSSM parameter space which 
provides the best fit || to the electroweak precision data, all four neutralinos are light enough 
[while the charginos are heavier than the present experimental bounds || from LEP1.5] for 
the top quark to decay into the lightest stop squark plus one of the four neutralino states. 

In the minimal supersymmetric extension of the Standard Model (MSSM), the decay 
t — ► tX\ could have a sizeable branching fraction?] when it is allowed kinematically [[J . In a 
large area of the SUSY parameter space, this decay can have a signature similar to that of 
the standard decay t — > bW + — > bl + u or 6+2jets M. Indeed, the stop could mainly decay 
into the lightest chargino and a bottom quark, and the chargino will subsequently decay 
into a neutralino and a virtual W boson which then goes into l + v or 2— jets. This gives the 
same topology as the t — > bW + channel, except for the large amount of missing energy in 
the SUSY decay mode due to the escaping neutralinos, which leads to a softer transverse 
momentum spectrum of the charged leptons. Such decays could have an impact on the top 
quark studies at the Tevatron ||, altering for instance the reported production cross section 
values^ since these events would not all pass the experimental cuts designed for the standard 
t — > bW + decay. 

A detailed study of the t — > txi decay mode is therefore mandatory. To have a more 
reliable prediction for the decay branching ratio, the QCD corrections have to be taken into 
account. The one-loop QCD corrections to the two other decay modes of the top quark 
in supersymmetric theories, t — > bW + and t — > H + b, can be found in Refs.|| and || 
respectively. In this article, we present the 0(a s ) QCD corrections to the top decay into a 
stop squark and a neutralino. These results can be easily generalized to the possible decay 
into a light sbottom squark and a chargino state. 

-"^The present agreement between the top quark production cross section at the Tevatron 0] and the 
Standard Model prediction still allows appreciable t — > t\i decay rates. A detailed analysis of the area of 
the parameter space where these decays are allowed or disfavored by Tevatron data is beyond the scope of 
the present paper. 

2 For such light stop squarks, the pair production pp — > it* is possible at the Tevatron and will also 
contribute to the event samples; for a recent discussion of this possibility see for instance Ref.Q. 
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2. Born Approximation 

In the MSSM, the decay width T(t — > Ux®) [where i — 1, 2 are for the stop and j = 1-4 for 
the neutralino states] will depend on five parameters: the left- and right-handed scalar stop 
masses [which in general are taken to be equal], the Higgs-higgsino mass parameter /i, the 
soft-SUSY breaking trilinear coupling A t , the ratio of the vacuum expectation values of the 
two Higgs doublet MSSM fields tan f3 [which fix the stop masses and the mixing angle] as 
well as the wino mass parameter M 2 [which together with \x and tan (3 fixes the neutralino 
masses]; the gluino mass is related to the parameter M2,rrig ~ 3.5M2, when the gaugino 
masses and the three coupling constants of SU(3)xSU(2)xU(l) are unified at the Grand 
Unification scale. 

In the Born approximation, the partial widths for the decays t — * tiXj [we will drop for 
simplicity the indices for the neutralino masses, m x = m x o\ are given by 



a 



a\ a l R (Am t m x e x ) + (a l L 2 + a l R 2 ) (m 2 t -m\ + m x ) \ A 1/2 (m 2 , m 2 , m\ 



'V, 



where X(x, y,z) = x 2 + y 2 + z 2 — 2 (xy + xz + yz) is the usual two-body phase space function 
and e x is the sign of the eigenvalue of the neutralino xi t he couplings a l L R which form the 
Born xtti vertex 



le 



with Pl,r the chirality projectors Pl,r 



bm t 




+ c L 




a l L P L + a R P R 
(1 =F7 5 )/2, are given by 

= bm t 



'R 



'R 



(2) 




\/2M w sin (3sw 



Na 



cl 
cr 



V2 
-V2 



3 " 3 
2 

3 



s 



w 



1 

2 3 
N' 

3c w 32 



Cy/Sw 



(3) 



Here 9 is the stop mixing angle [which can be expressed in terms of the Higgs-higgsino SUSY 
mass parameter /i, tan f3 and the soft-SUSY breaking trilinear coupling A t ] with se = sin 6, 
eg = cos 9 etc.; s^y = 1 — c 2 ^ = sin 2 ^^. is the diagonalizing matrix for the neutralino 



states 10 and 



N' n = c w N n + s w N j2 , N' j2 = -swNj! + c w N j2 



(4) 



The masses of the two stop squarks and their mixing angle are determined by diagonalizing 
the following mass matrix: 
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where M LR in the off-diagonal term reads M LR — A t — // cot (3. 

Assuming that the MSSM charged Higgs boson, H + , is heavier than the top quark 
so that the decay channel t — > H + b is kinematically closed, the branching ratios for the 
t — > hx^ decay modes [t 2 is the lightest stop squark in the notation introduced above] are 
shown in Fig. 1 for tan/? = 1.6, a value favored by theories with b-r Yukawa coupling 
unification. The top quark mass is fixed to m t = 180 GeV. In Fig. la, we fix the wino mass 
to M 2 = /i = 150, 200 and 250 GeV and plot the branching ratio for the decay into the stop 
and the lightest neutralinos as a function of the t 2 mass [obtained by varying the parameter 
At]. In Fig. lb, M 2 is kept fixed at M 2 = 400 GeV, and two values have been chosen for 
fi, —100 and 150 GeV. Whenever the decay in more than one of the neutralinos is possible, 
the branching ratio has been summed. As can be seen, it can reach the level of several ten 
percent, and can even dominate over the standard t — > bW + decay mode. 

In Fig.lc, we choose the parameter /i to be negative and equal to the wino mass M 2 = 
— fi = 50 GeV; this area of the MSSM parameter space, together with a rather light right- 
handed stop squark, mi < 100 GeV, provides the best MSSM fit to the high-precision 
electroweak data, and enhances the ratio Rb = T(Z — > bb)/T(Z hadrons) to the level 
where it is almost consistent with its experimental value 0. The branching ratios of the 
four decays t — » t + x°j with j =1-4 [the masses of the neutralinos can be read off from the 
thresholds, and the lightest chargino in this scenario has a mass m x + ~ 70 GeV] are shown 
together with their sum. They are also large, and can dominate over the t — > bW + branching 
ratio. 



3. QCD corrections 

The QCD corrections to the decay width, eq.(l), consist of virtual corrections Figs.2a-d, and 
real corrections with an additional gluon emitted off the initial or final state, Fig.2e. The 
0(a s ) virtual contributions can be split into gluon and gluino exchange in the t—t—x vertex, 
and the top and stop wave function renormalizations. The renormalization of the t—t—x 
coupling is achieved by renormalizing the top quark mass and the stop mixing angle. We will 
use the dimensional reduction scheme Jll| which preserves supersymmetryQ to regularize the 



ultraviolet divergencies and we introduce a fictitious gluon mass A to regularize the infrared 
divergencies. The sum of the virtual corrections obtained by taking the interference between 
the Born and QCD corrected amplitudes, and the real correction, obtained by squaring the 
sum of the two diagrams in Fig.2e, is finite as it should be. 

3.1 Virtual Corrections 

The QCD corrections to the t-ti~x vertex, eq.(2), can be written as 

*T = fe-i £ [gi l p l + gi r p r ] (6) 

j=g,g,mix,ct 



3 Because the \it coupling is renormalized by QCD, the dimensional regularization and dimensional re- 
duction schemes give slightly different results. 
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where G g , G~, G l mix and G l ct denote the gluon and gluino exchanges in the vertex, and the 
mixing and counterterm contributions, respectively. 

The gluon exchange contributions [Fig. 2a] can be written as 

G l,L = a L F l + a R F 2 ) G l,R = a R F l + a L F 2 ( 7 ) 

with the form factors F[ 2 given by 

F{ = B (ml,X,m i J+2m 2 t C -2ml(C 11 -C 12 )+2m 2 x C 11 

F\ = -2m t m x (C + C 11 ) (8) 

where the two and three-point Passarino-Veltman functions, B and C = C (m 2 ,m~ ,m 2 , 
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2 A 2 , m~ ) can be found in Ref . [[12] . 



The gluino exchange contributions [Fig2b], are given by 
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(summed over) and i the outgoing squarks], the form factors F lk are given by 
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The mixing contribution is due to the diagrams of Fig. 2c where the top-gluino exchange 
and the quartic scalar interaction contributions switch the Born ytt\ vertex to the 0(a s ) 
Xtt 2 one and vice versa; it is given by 
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3.2 Counterterms 

The counterterm contributions in eq.(6) are due to the top and stop wave function renor- 
malizations [Fig. 2d] as well as the renormalization of the top quark mass m t and the mixing 
angle 6, which appear in the Born coupling, 

G]£ L = ^a]; 2 (SZl + SZ^ +b{s e ,c e }5m t + bm t {cg,-sg}56 - c L {s e ,c e }56 

G l £ R = -a R ^ (5Z R + SZ il2 ^ +b{c g ,-s e }8m t -bm t {se,ce}S9 + c R {c e ,-sg}Se (12) 

In the on-shell scheme, the quark and squark masses are defined as the poles of the 
propagators and the wave-function renormalization constants follow from the residues at 
the poles; the corresponding counterterms are given by 
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where the self-energies £ and their derivatives up to a factor a s /3ir which has been 
factorized out, are given in the dimensional reduction scheme by 
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Finally, the mixing angle counterterm is chosenF] in such a way that it cancels exactly the 
mixing contributions G % mix of eq.(ll) [following Ref.fTS]]], 

Am t m g c 2 eB Q (m\,m u m g ) + c 2e s 2e (A (m 2 2 ) - A {m 2 J) (15) 



59 



2 2 

m~ — m~ 

*1 *2 



4 Had we chosen the MS scheme, i.e. subtracting only the poles and the related constants in eq.(10), we 
would have been left with contributions which increase linearly with the gluino mass. 
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In practice, this reduces to simply discarding the contributions of the diagrams Fig.2c. 
The complete virtual corrections to the t -—>■ tiX decay width will be then given by 



a a s 
3mf 47r 



Re [K - m| + m 2 )K G\ + a R 
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These corrections are infrared divergent; the divergencies are cancelled after adding the real 
corrections. 



3.3 Real Corrections 

These corrections are obtained by adding the amplitudes of the diagrams Fig.2e with real 
gluon emission and squaring the resulting amplitude; the corrected decay width is then 
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4. Results and Discussions 

The 0(a s ) partial widths of the decays t —>■ t 2 x normalized to their Born expressions, are 
shown in Figs. 3 as a function of the lightest stop mass for the scenarios discussed in section 
2. For the numerical analysis, we choose mt = 180 GeV for the top quark mass and for 
the strong coupling constant we take the value a s = 0.118. In Fig. 3a, the corrections are 
shown for M 2 = fi = 150 (solid), 200 (dashed) and 250 GeV (dotted); they are negative 
for relatively small values of the gluino mass [solid line, m g ~ 530 GeV] and positive for 
larger values of m g [dotted line, m g ~ 900 GeV]; for intermediate values of the gluino mass 
[rrig ~ 700 GeV, dashed line], they vary between -1% and +1%. In Fig. 3b M 2 is again fixed 
at 400 GeV and fi = —100 (solid) and 150 GeV (dashed). The QCD corrections are small 
in general, but they can reach the level of 15 % for some stop masses. 

The QCD corrections to the top decays into the lightest stop squark and one of the four 
neutralinos, t — > tx®, in the scenario where M 2 = —fx = 50 GeV are shown in Fig.3c. Here 
again, the corrections are rather small being of the order of 2 — 12 % in all cases. As in 
Figs.la-c, the sum was computed if more than one neutralino contributes. 

For M 2 = 400 GeV and hence for large gluino mass, m g ~ 1.4 TeV, the QCD corrections 
are negative and increase in absolute value, reaching the level of 12% for light stop squarks. 
This means that the gluino does not decouple for large m g . There is a logarithmic dependence 
of the correction on the gluino mass which is visualized in Fig. 4, where the correction is 
plotted against the gluino mass for ft = —100 GeV and mi = 50 GeV. This logarithmic 
behaviour is due to the wave function renormalization and is a consequence of the breakdown 
of supersymmetry as discussed in Ref.]15|], where a similiar feature appears in the decay of 
a squark into a massless quark plus a photino. 

5. Conclusion 

We have calculated the 0(a s ) QCD corrections to the top decay into the lightest stop 
and a neutralino in the minimal supersymmetric extension of the Standard Model. These 
corrections can be either positive or negative and increase logarithmically with the gluino 
mass. For gluino masses below 1 TeV, they are at most of the order of ten percent and 
therefore, well under control. 

In the case where the top quark is lighter than its scalar partners, the reverse process 
t —>■ txo might occur. This decay, together with the decay t — > bx + , would be the dominant 
decay mode of stop squarks if they are lighter than m t + rrig. The QCD corrections to these 
processes are under investigation [ TO . 
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Figure Captions 

Fig. 1: Branching ratios for top decays into the lightest stop and neutralinos as a function 
of the stop mass in the Born approximation; tan/3 is fixed to tan/? = 1.6 and the 
top quark mass to m t = 180 GeV. (a): M 2 — fM — 150 (solid line), 200 (dashed line) 
and 250 GeV (dotted line); (b): M 2 = 400 GeV and // = -100 GeV (solid line) and 
M 2 = 400 GeV and [i = 150 GeV (dashed line); (c): M 2 = -fi = 50 GeV the thin 
solid, dashed, dotted and dashed-dotted lines correspond respectively to the branching 
ratios of t — > t 2 x° with j = 1-4 (the neutralinos are ordered with increasing mass), the 
thick solid line corresponds to the sum over all contributing neutralinos. 

Fig. 2 : Feynman diagrams relevant for the 0(a s ) QCD corrections to the decay of the top 
quark into a stop squark and a neutralino. 

Fig. 3: Relative size of the 0(a s ) QCD corrections to the decay rate t — > t\ as a function of 
the lightest stop squark mass. We use the same inputs as in Figs.l. 

Fig. 4: The size of the QCD corrections as a function of the gluino mass rrtg with M 2 = 400 
GeV, /x = -100 GeV and m h ~ 50 GeV. 
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